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Abstract 

We consider a Bradley-Terry model in random environment where each 
player faces each other once. More precisely the strengths of the players 
are assumed to be random and we study the influence of their distribu¬ 
tions on the asymptotic number of potential winners. First we prove that 
under mild assumptions, mainly on their moments, if the strengths are 
unbounded, the asymptotic probability that the best player wins is 1. We 
also exhibit a sufficient convexity condition to obtain the same result when 
the strengths are bounded. When this last condition fails, the number of 
potential winners grows at a rate depending on the tail of the distribu¬ 
tion of strengths. We also study the minimal strength required for an 
additional player to win in this last case. 

Keywords and phrases : Bradley-Terry model, random environment, paired com¬ 
parison. 

AMS 2010 subject classification : 60K37, 60G70, 60K40 


1 Introduction and Results 

We consider here a model of paired comparisons which may be used as a proxy 
for sport competition, chess tournament or comparisons of medical treatments. 
A set of N players (teams, treatments, ...) called {1,..., A } face each other 
once by pairs with independent outcomes. When i faces j, the result is described 
by a Bernoulli random variable that is equal to 1 when i beats j and 0 
if j beats i (hence Xij = 1 — Xj^i). The final result is given by the score 
Si = '^j^i Xij of each player that is its number of victories. We call winner 
every player that ends up with the highest score. 

To each player i is assigned a positive random variable R modeling its intrinsic 
value, that is its ’’strength” or its ’’merit”. Given Yf = (R,..., Rv), the 
distribution of {Xij)i<i^j<N follows the Bradley-Terry model: all matches are 
independent and 

VI < * < J < A, P (W., = llVf) = . (1) 
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The distribution of Yi is chosen as follows. Let = {Ui,... ,Un) denote 
i.i.d. random variables; to avoid trivial issues, suppose that all Ui are almost 
surely positive. For any i S { 1,..., A^}, Vi denotes the i-th order statistic of 
the vector U^: the larger the index of a player is, the ’’stronger” he is. 
Bradley-Terry model has been introduced independently by Zermelo [25j and 
Bradley and Terry [3], it was later generalized to allow ties [siiin] or to incor¬ 
porate within-pair order effects [7] , see [2] for a review. Despite its simplicity, it 
has been widely used in applications for example to model sport tournaments, 
reliability problems, ranking scientific journals (see [5] or more recently [1] for 
references). Bradley-Terry model has also been studied in statistical literature, 
see for example and references therein. 

Nevertheless, Bradley-Terry model has rarely been associated to random envi¬ 
ronment models (see however |22j l and, to the best of our knowledge, never from 
a strictly mathematical point of view. The addition of a random environment 
seems however natural as it allows to manage the heterogeneity of strengths of 
players globally, without having to look at each one specifically. It is a method 
already used fruitfully in other areas such as continuous or discrete random 
walks (see [24] or [11] for recent presentations). Our problem here is to under¬ 
stand how the choice of the distribution for the strengths of players influences 
the ranking of the players. In particular, does a player with the highest strength 
ends up with the highest score? And if not, what proportion of players might 
win? These problematics are related to the articles detailed below. 

• Ben-Naim and Hengartner [T] study the number of players which can win a 
competition. These authors consider a simple model where the probability 
of upset p < i is independent of the strength of players: 

VI < i j < Af, P(Nij = 1) =pli<j-I-(1-p) lj>j . 

For this model, they heuristically show with scaling techniques coming 
from polymer physics (see 0) that, for large N, the number of potential 
champions behaves as y/N. In the Bradley-Terry model in random envi¬ 
ronment, Theorem [3] shows that the class of possible behaviors for this set 
is much richer. 

• Simons and Yao [21] estimate the merits based on the observations 
of (A'ij)i<i<j<Ar. They prove consistency and asymptotic normality for 
the maximum likelihood estimator. It is interesting to notice that this 
estimator sorts the players in the same order as the scores Si (see [12]). 
In particular, the final winner is always the one with maximal estimated 
strength. Theorem [T] shows that usually this player has also maximal 
strength when the merits are unbounded but it is not always true in the 
bounded case, see Theorem [2l 

Throughout the article, U denotes a copy of Ui independent of U)^, Q de¬ 
notes the tail distribution function and supp Q its support, P denotes the an¬ 
nealed probability of an event with respect to the randomness of Y^ and 
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while Py denotes the quenched probability measure given Yf, 
that is P ( • \Yi ). In particular, 

VI < i < J < iV, Py (X,,, = 1) = . 

We are interested in the asymptotic probability that the ’’best” player wins, that 
is that the player N with the largest strength Vn ends up with the best score. 
The following annealed result gives conditions under which this probability is 
asymptotically 1 when the number of players TV —> oo. 

Theorem 1 Assume that there exist /? G (0,1/2) and xq > Q in the interior of 
supp Q such that is convex on [ccojOo) and that E [f/^] < oo. Then, 

P ( the player N wins) > P f Sn > max Si ) - > 1 . 

l<i<7V-l J N^oo 

Remark 1 When the support of the distribution ofU is K.+ , the convexity con¬ 
dition is not very restrictive as it is satisfied by standard continuous distributions 
with tails function Q{x) ~ e~^ , Q{x) ~ x~^ or Q{x) ~ (loga;)"'”. The moment 
condition E [[/^] < oo is more restrictive in this context but still allows for 
natural distributions of the merits as exponential, exponential of Gaussian or 
positive parts of Gaussian ones. It provides control of the explosion of maximal 
strengths. It is likely that it can be improved, hut it is a technical convenience 
allowing to avoid a lot of tedious computations. 

When supp Q is finite, we can always assume by homogeneity that it is included 
in [0,1], since the distribution of given Yi is not modified if all 

Vi are multiplied by the same real number A. The moment condition is always 
satisfied and the only condition is the convexity one. This last condition forbids 
an accumulation of good players with strength close to 1. 

Let us investigate now the necessity of the convexity condition. For this purpose, 
suppose that Q{1 — u) ^ u°‘ when it —?► 0, then the convexity condition holds 
iff a > 2. To check the tightness of the bound 2, we introduce the following 
condition. 

Assumption 1 The maximum o/suppQ is 1 and there exists a € [0,2) such 
that, 

log(3(l — It) = alog(it) + o(logu) whenu^O. (A) 

Let us stress here that Q may satisfy m even if it is not continuous. Moreover, 
a is allowed to be equal to 0, in particular, (3(1) may be positive. Notice that 
some standard distributions satisfy Assumption for example the uniform 
distribution satisfies with a = 1, the Arcsine distribution satisfies it with 
0=1/2 and any Beta distribution B(a, b) satisfies it as long as the parameter 
b < 2 with a = b. 

The next quenched result studies, under Assumption the size of the set of 
possible winners. 
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Theorem 2 For any r G K+, let Gr = { |"A^ — r] + 1,..., iV} denote the set of 
the [rj best players. If holds, for any 7 < 1 — a/2 then, P almost-surely. 


Fv {none of the N'^ best players wins) = Py ( max Si < max S'i I —?► 1 . 

yiGGAT'Y J 

For any 7 > 1 — a/2 then, P almost-surely, 

Py ( one of the N"^ best players wins ) = Py ( max Si < max S'i I —>■ 1 . 

yi^G^'y i^Gm J 

Remark 2 The first part of the theorem shows that, when Q{1 — u) ~ u°‘, 
with a < 2, then none of ’’best” players, for any 7 G ( 0,1 — a/ 2 ) wins the 
competition. In particular, the ’’best” one does not either. In this sense, the 
bound 2 in the asymptotic development of Q around 1 is tight. 

The second result in Theorem H shows the sharpness of the bound 1 — a /2 in 
the first result. Heuristically, this theorem shows that under Assumption 
“/2 players can be champion. 

Under Assumption the best player does not win the championship. There¬ 
fore, we may wonder what strength v^v+i an additional tagged player iV -|- 1 
should have to win the competition against players distributed according to Q. 
The following quenched result discusses the asymptotic probability that player 
A" -|- 1 wins depending on its strength utv+i- To maintain consistency with the 
previous results, we still use the notations 


N 


S,= Y, fori € {l,...,iV + l} 




With this convention, Sn+i describes the score of the player -|- 1 and the 
score of each player i G {1,...,A^} is equal to Si -\- W,Ar-i-i- 

Theorem 3 Assume and let 



and 



//lim infAr_,.oo 


t'W + l —1 

(-N 


> 1, then, F-almost surely 



//lim supjv_^oo ^ < 1, then, F-almost surely 
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Remark 3 This result shows a cut-off phenomenon around 1 + eAr for the 
asymptotic probability that player + 1 wins. 

It is interesting to notice that, for a given a, ejsi is a non increasing function of 
Du- Therefore, when U is stochastically dominated by U', that is P(C/ > a) < 
P([/' > a) for any a € [0,1], we have djj < ■djj', hence . In other 

words, it is easier for the tagged player to win against opponents distributed as 
U' than as U even if the latter has a weaker mean than the former. This result 
may seem counter-intuitive at first sight. In the following example in particular, 
it is easier for the additional player to win the competition in case 1 than in 
case 2, since both distributions satisfy © with a = 0 . 

1. All players in {1,..., N} have strength 1. 

2. The players in {1,..., N} have strength 1 with probability 1/2 and strength 
1/2 with probability 1 / 2 . 

Actually the score of the tagged player is smaller when he faces stronger oppo¬ 
nents as expected, but so is the best score of the other good players. 

Remark that the first theorem is an annealed result while the others are quenched. 
Indeed, the first theorem requires to control precisely the difference of strengths 
between the best player and the others when all the players are identically dis¬ 
tributed, this seems complicated in the quenched case. This problem does not 
appear in the other results: for example, in Theorem [31 the strength of the 
tagged player is deterministic and the strengths of others are bounded by 1 . 

The remaining of the paper presents the proofs of the main results. Section [2] 
gives the proof of Theorem [T] and Section [3] the one of Theorems [2] and [31 

2 Proof of Theorem [1] 

Denote by Zjs; = maxjg r ^ } >5'i- The key to our approach is to build 

random bounds s^ and depending only on such that, 

^{Sn>s^)^1, ^{Zn<z^)^1 and ¥{s^>z^)^l. ( 2 ) 

It follows that, 

P(S'iv > Zn)>¥{Sn > s^, Zn < z^, s^ > z^) 

>1-V{Sn <s^)-F{Zn> z^)-F{s^ <z^) ^1 . 

The construction of s^ and z^ is the subject of the next subsection, it is 
obtained thanks to concentration inequalities. The concentration of Sn is easy, 
the tricky part is to build z^. First, we use the bounded difference inequality 
to concentrate Z^ around its expectation. The upper bound on its expectation 
is given by the sum of the expected score of player N — 1 and a deviation term 
that is controlled based on an argument used by Pisier [19]. Finally, the control 
of P (s^ > derives from an analysis of the asymptotics of Vn-i and Vn- 
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2.1 Construction of and 

The expectation of the score S'at of the best player is given by 


N-l 


Ev[Sn]=J2 


Vn 


i=l 


Vn + V, 


and the concentration of Sjv is given by Hoeffding’s inequality, see |15j : 

/ N-l 


Py Sat < ^2 


Vn 


Vn + Vi 


< e 


Hence, the first part of (I2|) holds for any un —>■ oo with 


Af-l 




= E 

i=l 


Vn 


Vn + V, 


- 


UN ■ 


(3) 


The following lemma implies the concentration oi Zn around its expectation. 
As we will use it in a different context in other proofs, we give a slightly more 
general result. 

Lemma 4 Let I C [A^] and let Z = max^g/ Si. For any u > 0, 


Py ( Z >¥.v [Z] + y —u I < e 


N 


and 


Py 


{^Z<¥.v[Z] 



Proof: The proof is based on the bounded difference inequality recalled in 
Theorem [T^ of the appendix (see [3], [H] or m)- To apply this result, we have 
to decompose properly the set of independent random variables {Xij)i<i^j<N- 
To do so, we use the round-robin algorithm which we briefly recall. 

First, suppose N even. Denote by tr the permutation on {1,..., N} such that 
cr(l) = 1, (7{N) = 2 and (j{i) = i + 1^ it 1 < i < N and define the application 


A:{l,...,iV-l}x{l,...,iV}^{l,...,iV} 

by A{k, i) = {^N -I- 1 — (i)). Then, for any k G {1, ..., — 1}, 

A(k, ■) is an involution with no fixed point and for any i G {1,..., A^}, A{-, i) 
is a bijection from {1, ... ,N — 1} to {1, ..., A^} \ {*}. The first variable of 
function A has to be thought as “steps” in the tournament. At each step, every 
competitor plays exactly one match and A{k, i) represents here the opponent 
of player i during the step. We denote by Z^ the variables describing the 
results of the step, that is Z^ = i < A{k,i))- The variable Z 

can be expressed as a (measurable) function of the Z = ^(Z ^,..., Z^~^). 
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Moreover, for any k = 1,..., N — 1 and any ..., ^,z^ in { 0,1 the 

differences 


|^'(2 


^N-1 




^N-l\ 


are bounded by 1. If N is odd, we only have to add a ghost player and Z can be 
expressed in the same way as a measurable function of N independent random 
variables with differences bounded by 1. Therefore, in both cases, the bounded 
difference inequality applies and gives the result. □ 

It remains to compare the expectations of and of Sn- This requires to 
control the sizes of Vv-i and Vn. Recall that Vi is the order statistics of 
the vector = (C/i,..., U_n), that is 

Vi = min {Uk I 3/ C {1,. .., IV}, |/| = i and VI € I,Ui < Uk] ■ 

ke{l,...,N} '' ' ■’ 


Then the sets { Vi,..., Vn } and {Ui,... ,Un}, counted with multiplicity, are 
equal which guarantees that, for any function /, /(^) = f{Ui)- 

Let Q~^ denote the generalized inverse of Q: for y G (0,1), 

= inf {a; G K+, Q(a;) < y} . 

Remark that the convexity assumption implies that, if M is the supremum of 
suppQ, the function Q is a continuous bijection from [xo,M) to (0,(3(xo)] such 
that Q{x) = 0 so, on (0, Q{xo)], Q~^ is the true inverse of Q. 

Lemma 5 For any function h defined on K+ such that lim_|_oo h = +oo, let 

0 otherwise 

= {a% < Vn- 1 <VN<b%} . 

Then lim P ( A %) = 1. 

Proof: If h{N)/N > 1, P {Vn-i < a,%) = 0 so the lower bound is trivial. If 
h{N)/N < 1 and h{N) > I, since x A Q{xo) < QiQ~^{x)) < x, 

P ( Vat-i < = (1 - Qia%)f + iV(l - Qia%)f-^Qia%) 


Hence, P (Vat-i < «jv) 0. Moreover, for any N such that l/{Nh{N)) < xq, 

/ 1 \ ^ 

P(lV>6'V) = l-P(C/<6^)^ = I-(^I-^^^j ^0. 

□ 
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Lemma 6 There exists a non-increasing deterministic function y —^ "niu) oti 
R+ such that lim+oo rj = 0 and limjv-^oo P [Bn ) = I, where = | < r]{N) 


Proof: Since E[U^] = /“ yQ{y)dy < oo, lim^j-^+oo /“ yQ{y)dy = 0. As Q is 
non increasing, 



ydyQix) < 



yQ{y)dy -)> 0 . 


Therefore Q{x) = o(l/a:^) when x —>■ +oo. This implies that Q~^{y) = o(l/y^) 
when ?/ —>■ 0 and there is a non-decreasing function y u{y) defined on K+ 
such that limo u = 0 and Q~^{y) < u{y)/y/y. For any N large enough, choosing 

y ~ N _ 

Q-' < VN^a/N)- 

Setting r]{y) = 'Juiftjy')^ Lemma [5] used with h{x) = 1/ ^yu{l/x) gives the 
result. □ 


We also need the following result: 


Lemma 7 Define En{V) 


N ^ Vn-1 + F, 
2=1 


and 


Cn = \ [[/] < En{V) < 2E [U] 


Then, limAr_>oo E ( Fat ) = 1. 

Proof: Remark that for i G {1,...,A^—1}, 1 < Id- ^ < 2, hence 


N-l 


En{V) > 


2N 


i=l 




N 


i=l 


Vn 

m 


and 

JV-l N 

2=1 2=1 

Lemma ensures that V^/N converges in probability to 0 and therefore the 
proof is easily conclude using the Weak Law of Large Numbers. □ 

We are now in position to bound Ey [ Zn]- 


Lemma 8 For N large enough, on n Cn , 


Ev[Z n] < 


N-l 


E 


Fv-i 

Fv_i+F 



A^logA^ 

Fv-i 
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Proof: Write 


N-l 


Vn-1 


Z'^-Z^-Y, y ^ y - 
1—1 

For any A > 0, Jensen’s inequality and the argument of Pisier 

[Z'n] < j logEv 
1 

< ^ log ( 


give 


2=1 


\f C V^N-1 '^N-l 

= l Vjsr.i + V, 


Let S = where Xi are independent Bernoulli variables with respective 

parameters Vn-i/{Vn-i + Vi). Every Si is stochastically dominated by S, thus 


Ey [ Z'jsf ] < — log { NEv 


\ f O’ ^ — 1 - — 1 

V^ZTTV- 


Let 


\n = 


IVn-i logN 
2NE[U] 


for N large enough, on Bat fl Cn, Aat < 3/(8e^) < 1. Lemma [12] in Appendix 
evaluates the Laplace transform of Bernoulli distribution and gives here 


Ey 


An 


e 




Vn-1 ^ 1 

Vfj-l + Vi J 


< e 


^2 

■^N Z^i=l 


VN-lVj 

(VN-l+VP^ 


Therefore, 


□ 


Ey [^Af] — 


logfV 

Aat 


2XnN 

Vn-1 


E[C/] 



iVlogAf 

IV-i 


Lemmas |31ini|7| and |5| give the second part of (O for any uat —>■ oo with 


z 


N 


N-l 


E 


Vn-1 
Vn-1 + r, 




(4) 


To prove the third item of (I2|) , it remains to prove that there exists some un ^ 
oo such that the probability of the following event tends to one: 



Vn 

V^PVn 


Vn-1 \ 

r, + ) 


> 


/8E[17] 


iVlogA^ 

Vn-1 


+ 2 a/V 


Un 


(5) 
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2.1.1 Proof of ([5]) 

The proof relies on a precise estimate of the difference — V/v-i- 


Lemma 9 Let 


Dn 


Vn-1 < Vn 



Vn-1 

y/N 



Then P (Dn) ^ i as N ^ oo. 

Proof: In the proof, Cp denotes a deterministic function of /3 which value can 
change from line to line and F = 1 — Q denotes the c.d.f. oiU. 

Let r] be the function defined in Lemma [Bland denote h = rj~^. As Xq is in the 
interior of the support of Q, there exists a constant c < 1 such that Xq/c also 
lies in the interior of supp Q. Let a^, be defined as in Lemma [5] and consider 
the event 


Ttv = I (xo/c) V < Vn -1 < VN (1/2 A 'q{N)) A } 


According to Lemma [5] and Lemma IH limAr_>.oo P (-FV) = 1- Moreover, 


P{D%)<P{D‘i,f^FN)+P{F|,) 

= n^F,^{D%\VN-i)]+P{F^N) 


The cumulative distribution function of the random variable Vn given Vn-i is 
1 — Q/Q{Vn-i) and then 




Q(yN-i) 


On the event Fat, the convexity of Q gives 



and Vn/Vn is smaller than 1/2, thus 


P{D%\Vn-i)<Cp 


V^-Sf'{Vn-i) 


( 6 ) 


m^-p)/^Q{VN-i) 
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Moreover, by convexity of the function F 'is non-increasing, 

hence, by Holder’s inequality, for any x > a;o/c. 


^2-2/3 


Q(x)l/2+/3 


Q(7/)1/2+/3 


dy 






< Cp log 


Q{cx) 

Q{x) 


1/2-K/3 


As seen in the proof of Lemma [U limy_>.oo y^Q{y) = 0, hence Q{cx) < Cp^xjx^ 
for some constant C^^\. Therefore for any a; > xo/c. 


< Cp ( log 


Q(x)l/2+/3 


g/3,1 

3^Q(x) 


1/2+13 


The function g{x) = (log j is upper bounded, this 

yields the following inequality: 

x^-^^/^F\x) < CpQ{xf/^+^P/^ . 

This last bound applied to x = Vm-i combined with (|6l) leads to: 

¥{D%\Vn-i)< - C;3(14f-i)^/ 

iV(l-«/2(g(y^_l))l/2-3/3/4 
Now on Tjv the following bounds also hold: 


Vat-i < VNr]{N) and g(VAr_i) > 


y{N)^ 

N 


hence, 


lF:,lPiD%\VN-i)<CpyiNf^"/^ . 

We conclude the proof by integration of the last inequality. □ 

Let h(x) = x/2 and consider the event 

Gat = n Bjsf n G^r n Dj^ 

= {vA-. > <3-‘(i/2). ^ - V 7 ._. > v„ ( 1 ^)' 

^^<£;jv(' 1 /)< 2 E[G]|. 
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According to Lemmas [5l IHl 0 and O P (Gat ) converges to 1 when N ^ oo so 
we only have to prove that (l5|) holds on the event Gn- And on Gat, 


N-l 

E 

i=l 


Vn 


Vn-i 


Vi + Vn Vi + VAf-i 


^ Vn — Vn- 1 Vi 

^ Vn- 1 + Vi 
2 = 1 


2Vn 


> 


^y ^ ^ E[y]^/yv y 

Vn-1 J 2 - 8 [ Vn-1 ) 


Thus, for N large enough, on Gat, 

Y^Af—1 / Vpf _ Vw-1 \ 

l^i=l Vi+Viv Vi+VN-I ) 


Vi+VN Vi+VN-1 

8 E[G]Sf 


> 


> 




83 


(Viv-i) 


VioglV 


"I'^I(< 3 -( 1 / 2 ))''« 4 S ,>2 (T) 


83 


n/IoglV 


and 

iV-l 

E 


Vn 


Vn-1 


14 + Vn V + Vn- 


A > f A^y > ^ 

1 J 8 y Vat-i J 8 r]{N)P 


Hence, for a constant c small enough and un = c/ {rj^N) on Gat, 


N-l 


2\J Nun < 77 E 


i=l 


Vjv 


Vn-1 


+ Vn + Vv-i 


( 8 ) 


Bounds d?]) and ([8]) imply ([5|) ; this concludes the proof of Theorem 0 


3 Proof of Theorem [2] and Theorem [3] 

Remark that in both theorems, each variable Si for z S { 1,..., A^} has the same 
definition, it corresponds to the score of a player with strength Vi playing against 
opponents with respective strength {Vj,j G A^} \ {i}}- Therefore in 

both proofs, the notation Zn = maxjgr atj Si represents the same quantity. 
We build bounds and depending only on such that 

Pv (s- < S’at+i < S^) ^ 1, Pv (zy < Zat < ^ 1, P-a.s. (9) 

and such that, when lim inf at-s-oo > 1, P-almost surely, for any N large 

enough, > 1 + z^, while when lim supAr_>.oo P-almost surely, 

for any N large enough, < z^. In the first case, it follows that, on 

{sy>i + zy}, 

Pv ( S'at+i >1-1- Zn ) > Py ( Sn+i > , Zn < 2+ ) 

> 1 - Py {Sn+1 <s^)-¥v{Zn> z^) . 
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The result in the second case is obtained with a similar argument. This will 
establish Theorem 

For Theorem [21 given 70 < 1 — a/2 < 71 , we build random bounds and Zi 
depending only on such that P-almost surely, 


Pv 


( max 
y 260^70 


Si < z. 


N 




P 


( max Si < z^ 


-> 1 . 


and 

P (liminf { } ) = P (liminf { zf < } ) = 1 . 

On {z^ <z5}, 


( 10 ) 


i^G 


Pv ( max Si < Zm ) > Pv ( max Si < Zq , z^ < Zm 




> 1 — Py ( max St> Zq] —Pv{z- < Zn) 
\ieGff~io ) 


On {zf < z(/}, 


Pv" ( max Si < Zjq ) > Py ( max Si < zf^, z^ < Zn | 

Vi^Gjv7i / \*^G,,7 i / 

> 1 — Py 1 max Si > zf^ ) — Py (z(/ < Zn ) . 
V*^G„71 V 

Together, these inequalities yield directly Theorem |2l 

The construction of and will derive from the concentration of Sn+i given 
by Hoeffding’s inequality: for any u > 0, 


Py 


Py 


N 


Sn+1 < ^ 

^ i—1 

' N 

sn+1 > 


i=l 


VN+1 

Vi + vn+1 


VN+1 
Vi + uat+i 



( 11 ) 

( 12 ) 


We will now build the bounds Zq , z^, z^ and z//. To do so, we study the 
concentration of Z^, maxigCit^ Si and maxi^g^^ Si. The construction of these 
bounds is based on the same kind of arguments as the ones used in the previous 
section. The construction of z(/ requires a lower bound on Ey [Z^^] which is 
obtained by comparison with the maximum of copies of the Si that are inde¬ 
pendent, see Lemma [TTJ 


3.1 Construction of , z^ and z^ 

LemmalUgives the concentration of Zjv, maxigG^.^p Si and maxj^( 3 ^.^j Si around 
their respective expectations which are evaluated in the following lemma. 
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Lemma 10 P almost-surely, 


Ey = NM 
Ey 


1 


l + i7 


max Si 

< fVE 

1 

.^ + u_ 




+ \/(2 — a)duN log N + o{y/W\ogN) 




where v = —— > 0. In addition, P almost-surely, 


Ey 


max Si 

< iVE 

1 

.1 + t^. 




yj2^QduN\ogN + o{y^N\og N) . 


Proof: 

Upper bounds. Define Z'n = Zff — i+Vk ■ iterated logarithm 

ensures that, P-almost surely, 


E 


1 

1 + 14 


E 


1 

1 + Uk 


fVE 


1 

1 + D 


' (^ViVlog 


N 


(13) 


To bound E [ Zn ], it is then sufficient to prove that 

Ey [Z'f^] < a/(2 - Q)'duN\ogN + o{\/N\ogN) . 

Let e > 0 and /|^ = {* s.t. 14 > 1 — }. By Jensen’s inequality and 

the argument of Pisier m, for any A > 0, 


Ev" [Z'n] < ^ log 


E + E 






Ey 



Let S = J2k=i where, given V(^, the Xk are independent Bernoulli variables 
with respective parameters 1/(1 + Vk) and S' = J2k=iZk, where the Yk are 
independent Bernoulli variables with respective parameters 

1 _ _/V“42+e 

1 _ N-^/2+e ■ 

The variable S represents the score obtained by a player with strength 1 playing 
against all the others, so it clearly dominates stochastically each Si. Likewise, 
S' represents the score obtained by a player with strength 1 — playing 

against all the others, so it dominates stochastically each Si, with * i ih- 
Therefore, 


Ey 


[Z'n] < 


(log 




„A(S-Ev[S]) 






■ (14) 


14 




































Let Xn = C ^/log N/N where C is a constant that will be defined later. For any 
1 < fc < 

1 + b-fc- y/N{l + VkY 


which gives following the upper bound: 

N 


Ev 




< 


k=l 


,XN{Yk-Ev[Yk]) 


—Xn 


Vk 

Yn (1+v^)^ ^ 


By Lemma [TSl for N large enough 


El 




N x2 

k=l 


_N±_ Vk 

Yn (i+Vfc)^ 




As the strong law of large numbers shows that, P-almost surely, 


N 


Uk 


we obtain 


Ev 


^ § (1 + 
a«(s'-ELit^) 


•v/log N 


/ 


\=o[ 

^ ) 


' — a.s. . 


(15) 


We turn now to the other term in the right hand side of (HI: using Lemma [T^ 
and the law of the iterated-logarithm, P-almost-surely, 


E 


'V 


„AK(S-ELiT4k) 


N 

REv 


,XN{Xk-E[Xk]) 


fc=l 


< e 


^^E^=iVar(Xfc)+Of ^ 


log^/^ N 


It remains to control |/^|. By d^, P-almost surely, 

P (c/ > 1 - = 7V-“/2-SeagO(logAr) ^ 

then it is easy to prove, applying Borel-Cantelli lemma, that 

|7^| = P - almost surely . 

Therefore, (fT^ . (fTOl) and (fT71) prove that, P-almost-surely 

Tu- rv/i/n /o , ^logiV , AfA at ^ , f\ogN 

[^n\ ^ (1 ~ H" cxc)— -1-- —iju H- o I —-- 

AN ^ \ an 


(16) 


(17) 
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Hence, choosing C = that is Xn = ^ we get 

¥.v[Z'j^]< \/(2 - a + 2ae)duN log N + o{^/N log N) P-a.s. . 


As the result holds for any e > 0 small enough, this gives the upper bound on 
Ev[ Zn ]. _ 


Proceeding as in the proofs of (fTdl) and (fTHI) , but choosing now Xn = 


7 q log N 
N'du ’ 


we get the upper bound for Ey [ maxigCjvTo ] ■ 

Applying (ITT)) with e = ly, we get that, P-almost surely, for N large enough, 


|GAr7i I = iVT-i = iVl-“/2+2a<^ > ^l-a/2+a^.gO(logAr) ^ |j^| _ 

Therefore, for any i ^ Gn~>-l , Vi < 1 — We can prove as in the other 

cases that P-almost surely. 


Ey 


N 


max 


s-E 


1 - l/iVl/ 2 -^ 
l-l/NV'^-’y + Vk 


= 0{VNlogN) . 


It remains to remark that, P-almost surely, by (USD and the strong law of large 
numbers. 


^ 1 - l/iVl/2- ^ 

1 - IINV^— + Vk ^ 1 + H 


N 




14 


^ (1 + Vk? 


= A^E 


1 


l + U 




This concludes the proof of the upper bound on Ey [maxj^g^^j^ 5^]. 

Lower bound on Ey [Zn]- Let us start with the following lemma which says 
that Zn stochastically dominates the maximum of independent copies of the 
variables Si. 


Lemma 11 For any a > 0 we have, ¥-almost surely, 

N 

Py < a) < J^Py (5'i < a) . 

i=l 

Proof: We proceed by induction, we provide a detailed proof of the hrst step, 
the other ones follow the same lines. Let X 2 ,i denote a copy of ^ 2 , 1 , independent 
of (Xij)i<i^j<N and let S 2 = ^ 2,1 + "^ 2,0 M 2 = Si\/ S 2 , M 2 = Si\/ S 2 , 

A = { maxi >3 S'i < a}. Write { Zn < a} = { M 2 < a}nA. Simple computations 
show that: 

Py ( Zn < a) =Py ^ I Zn ^ a j- ^ — Py ^ | M 2 = a j- D | M 2 = M 2 -I-11 D 

-f Py (^|m 2 = a-k i| n I A /2 = M 2 - 1 } n a) . 
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In addition, for x G {0,1}, 


{ N N 

Xi^2 = -^2,1 = X , X2,i > X + Xi^i 

i—3 i—3 

Now, recall that Xi ^2 and Xi ^2 are independent of X 2 ,i, Xi,i and A, 
then, for x G {0,1}, 



Pv" ^ I M 2 = a + a; I n I M 2 = M 2 + 1 — 2x | fl 

= Pv (^ 1,2 = -^ 2,1 = a;^ Py I i ^ X 2 ^i = a, ^ < a — a; I fl A 

V I i=3 i=a J 

As Py ^^ 1,2 = ^ 2,1 = 0^ = Py ^Ali _2 = ^ 2,1 = 1 ^ we obtain 
Py ( Zj^s < a) < Py ^ Zj^s < a ^ . 

□ 

Let = {i s.t. > 1 — } and S = X]^ 2 where the Xi are inde¬ 
pendent and Xi ^ Bil/ll + Vi/{1 — The variable S is stochastically 

dominated by any Si with i € It follows from Lemma [TT] that. 

N 

Py(ZAr<a)<]^Py(S', <a)< Py ( S', < a ) < Py ( S' < O ' . 

»=1 i6/^ 

For any e G (0,1 — a/2), denote 'Jn = \/{2 — a — 2e)'djjN log N. The previous 
inequality yields 

lEy [Zn] > {iN +py [S])Py {Zm — Ey [S] > 7iv) 

> (7,v+Ey [S])(I-Py(S-Ey [S] <7iv)'""') • 


Denote \n = y ■ By Lemma [121 for any u G M+ and any N such that 
uAtv ^ 1 ) 


log Ey 


p«An(S-E[S]) 


AT-l 

E 


i=l 

Y 


'N 




(1 + ^Y 


0{X%) 




t"-a.s. 


\ Vn J 

The last line is obtained thanks to the law of iterated logarithm. Hence, 


1 


lim 
Af-s-oo log N 


log Ey 


g ® y V log N 
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The same argument applied on the variables —Xi shows that the previous in¬ 
equality actually holds for any u S M. Therefore, using Theorem[T4]in Appendix 
with the sequence of random variables C,n = ^7^7 ^ ; 


lim inf : 


■ logPy (^ — Ey [S] > 77 V) > —1 + a/2 -I- e . 


N log TV 

In particular, since log |/^| ^ (1 — a/2) \ogN, for N large enough 
Py(5'-Ey[5'] <77 v)'^"' < 


< e 


-N 


i/4 


Since, by the law of iterated logarithm, 


1 


^l-l/^ + C/, l + Fi 
1 


> NE 


1 + U 


+ o{y/N log N) 


we obtain that, for any e > 0, 


Ey [^7v] ^ A^E 


1 

l + U 


+ yj{2 — a — 2e)'duN log N + o N log N ^ 


which concludes the proof. □ 


3.2 Conclusion of the proof of Theorem [3] 

Choosing u = log log N in (ED, a slight extension of the law of iterated logarithm 
gives that, P-almost surely, with Py-probability going to 1, 


N 


sn+1 > 


VN+1 


Vi + VN+l 


Nun 


= IVE 


VN+1 
U + VN+1 


+ o 


(^^/N\og 


N 


and 


N 

sn+1 < 

i=l 


VN+1 
Vi vn+1 



VN+1 

U -I- VA7+1 


-I- o 




Therefore, there exists —>■ 0 such that the first statement of (O holds P- 

almost surely with 


s7 = iVE 


U + Vn+1 


±e7^N\ogN . 
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By Lemma [3] with u = log log and Lemma ITUl P-almost surely, with Py- 
probability going to 1, 


Zpf = 7VE 


1 


= ATE 


U+1 

1 

U+1 


\/(2 — a)'duN log N + o ^ yjN log N ^ . 
such that the second statement of ([9|) ho 
+ -\/(2 - a)'&uNlogN + f?j^^jN\ogN , 


+ 

Therefore, there exists —>• 0 such that the second statement of (l9|) holds with 

+ 1 = A^E ^ 


+ \/(2 - a)'du N log N - e%\/NlogN . 


Moreover, 


E 


VN+1 

U + vn+1 


1 

1 + U 


{VN+I - 1)E 


U 

{U + l){U + VN+l) 


Hence, denoting the inequality + 1 is verified if 


(^^AT+l — 1)]E 


_ U _ 

(JJ + 1){U + fAT+l) 


> ^ \/(2 — a)'du + ^ 



that is if 


VN+1 


- 1 


E 


_ u _ 

(C/+l)((7+t;N + i) 


£N 


du 


> 1 




yj{2- a)'du 


(18) 


where cn = \/2 — a is the value appearing in the statement of 

Theorem [31 We now prove by contradiction that 


lim inf 
N—¥00 


Vn+1 - 1 ^ 

CAT 


_ U _ 

{U + l]{U+VM + l) 


■du 


> 1 . 


(19) 


Suppose there is a subsequence of wat+i (that we still call z;Ar+i) such that (fT3]) 
is not true. As liminf(?;Ar+i — l)/eAr > 1, it means that for N sufficiently large, 
t'Af+i > 1 + (5 for some (5 > 0. But in this case, the LHS of (IT3)) clearly goes to 
infinity as A^ —>■ oo. That contradicts our initial assumption and then (1181) is 
verified for N large enough. 

The proof that when liminf(?;Ar+i — l)/eAr < 1 follows the same 

arguments. 


3.3 Conclusion of the proof of Theorem [2] 

By Lemma ID with u = log log N and Lemma Uni P-almost surely, with Py- 
probability going to 1, 


max Si < A^E 


1 


i + u 


\/2'yo'0uN\ogN + o (^^/N log A^^ . 
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Since 70 < 1 — q:/2 , the first item of (fTUl) holds for N large enough with 


= NE 


1 

1 + U 



|)%lVlogiV . 


It is clear that Zq < z^ for N large enough since, by definition 70 + !—^ < 2 —a. 
By Lemma S] with u = log log N and Lemma HUl P-almost surely, with Py- 
probability going to 1 , 


max Si < NE 


1 

1 + U 




where ^ > 0 . 

enough with 


Hence, the second item of m holds for N large 


zf = A^E 


1 


1 + U 


which is clearly smaller than z^. 


Appendix 

To evaluate the various expectations of suprema of random variables, the follow¬ 
ing result is used repeatedly. Its proof is straightforward and therefore omitted. 


Lemma 12 Let X be a Bernoulli distribution with parameter p G [0,1] and 
a € [ 0 , 1 ], then 


1 + yP(l -P) < IE 


^a{X-ElX]) 


< e 


p(l-p)a^( i-|-i§^a) 


We also recall for reading convenience two well-known results. The first one is 
the bounded difference inequality (see Theorem 5.1 in [H]): 


Theorem 13 Let T" = Ti x • • • x Xn be some product measurable space and 
: T" —> K. fee some measurable functional satisfying the bounded difference 
condition: 

|4'(a;i,...,a:i,...,a:„) - 4'(a;i,...,yi,...,a;„)| < 1 

for all X G T”, y G T", i G {!,...,n}. Then the random variable Z = 
^'(Xi,..., X„) satisfies for any u > 0, 


P 


Z>E[Z] + 



< e 


and 


Z<E[Z] 



The second result is a simple consequence of Gartner-Ellis theorem (see Theorem 
2.3.6 in [To]) and of the Fenchel-Legendre transform of a centered Gaussian 
distribution. 
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Theorem 14 Consider a sequence ofr.v. (Cn)„gN ® deterministic sequence 
(a„)„gN oo such that for any m G R, 


lim — logE 

n—>-oo a^i ^ 


9 


2 


2 


Then, for any x > 0, 

1 

liminf — logP(C„ > x) > 

n—>-oo Qn 
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